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RECEIVER SYSTEM FOR MULTIPLE-TRANSMIT, MULTIPLE-RECEIVE (MTMR) 
WIRELESS COMMUNICATIONS SYSTEMS 

BACKGROUND OF THE INVENTION 

The present invention is directed to the deUvery of data via a wireless connection and, 
more particularly, to the accurate detection of data received via a wireless connection. 

In a wireless communications system, data is delivered from a transmitter to a receiver 
using radio waves or other electromagnetic waves. The waves, however, may be reflected off 
objects in the environment and scattered randomly while propagating from the transmitter to 
the receiver so that multiple images of the transmitted signal may arrive at the receiver via 
different scattered paths. The scattering, known as multipath scattering, can cause the images 
to arrive at the receiver at slightly different times and interfere destructively, thereby 
cancelling each other. Thus, multipath scattering may significantly impede the accurate 
transmission of data over the wireless connection. 

A known approach for achieving high spectral efficiencies is by transmitting the data 
using a base transceiver station (BTS) having multiple antennas and then receiving the 
transmitted data using a remote station having multiple receiving antennas. This approach is 
known as Multiple Transmit-Multiple Receive (MTMR). As an example, a single data 
stream is divided into multiple substreams, and an array of transmitter antennas is used to 





concurrently transmit the parallel substreams using the same frequency band. The transmitted 
substreams are then picked up by an array of receiver antennas. Each of the receiver antennas 
receives all of the transmitted substreams as a single, superimposed signal, rather than as a 
plurality of individual signals. A receiver located at the remote stations then demodulates the 
5 superimposed substreams to recover the original, transmitted substreams. 

The scattered images of the transmitted substreams, however, are also picked up by 
the array of receiver antennas. Because each of the substreams originates from a different 
transmitter antenna located at a respectively different point in space, each substream is 
scattered differently, i.e. independently, and is thus received independently at each receiver 
lQ,g antenna. When there is sufficient multipath scattering, the original substreams are not readily 
=p recovered. 

The scattered images of the transmitted substreams may arrive at the receiver antenna 
l^' at the same time, i.e. with no delay spread, known as flat fading or may arrive with at 
fy different times, i.e. with a delay spread, known as frequency selective fading. In the delay 
1 sen spread or frequency selective fading environment, currently transmitted symbols interfere 
with previously transmitted symbols, which is called inter-symbol intereference (ISI). 

A known signal processing algorithm, referred to as "Bell Laboratories Layered 
Space-Time" (BLAST), uses differences in the scattering to identify and recover the 
transmitted substreams. The substreams are first sorted according to the strength of each 
20 channel, and the strongest substream is detected and extracted from the total received signal. 
Then, the substreams are again sorted based on their signal strength, and the next strongest 
substream is detected and extracted from the total received signal. The sorting and extracting 
process is repeated until all of the substreams are determined. The BLAST algorithm 
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assumes that only flat fading is present and is therefore not suitable for delay spread 
envirorunents, such as in a mobile wireless environment. 

Thus, in the presence of ISI, the BLAST algorithm is not applicable. Moreover, while 
detecting the current strongest substream, the BLAST algorithm nulls the remaining 
5 substreams. This nulling is often imperfect because of channel estimation errors and, as a 
result, the performance of the receiver is degraded. 

It is therefore desirable to recover the received data streams in the presence of ISI and 
without degrading receiver performance. 

1 aa SUMMARY OF THE INVENTION 
^i; The present invention provides a wireless communication system for transmitting and 

Jh receiving a MTMR signal in which the superimposed substreams of the received MTMR 
s signal are recovered by estimating all of the substreams of the transmitted signal 
fU concurrently. 

1 In accordance with an aspect of the invention, transmitted data symbols estimated 

from a plurality of signal samples that are received by a plurality of receivers. The data 
symbols are part of a data stream that is divided into a plurality of substreams. Each of the 
receives signal samples includes one or more data symbols from each of the plurality of 
substreams. The plurality of signal s£imples are represented as a function of a plurality of 

20 transmitted data symbols and a plurality of channel response values. Each of the plurality of 
charmel response values represents a respective signal response for a transmission path of at 
least one of the transmitted data symbols. At least part of the plurality of channel response 
values includes channel responses caused by multipath scattering. A plurality of estimated 
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data symbols is defined as a function of the plurality of signal samples. The function satisfies 
a performance related criterion. A difference expression that represents a difference between 
a function of a plurality of estimated data symbols and a function of the plurality of signal 
samples is determined and is a sum of a plurality of terms. Values for each of a portion of a 
plurality of estimated data symbols are repeatedly selected such that one of the plurality of 
terms if minimized. Values for a further portion of the plurality of estimated data symbols 
are repeatedly selected such that a further one of the plurality of terms is minimized, where 
the further one of the plurality of terms is a function of the further portion of the plurality of 
estimated data symbols and the selected values of the plurality of estimated data symbols, 
until each of the terms is minimized. 

In accordance with this aspect of invention, the received signal samples each include 
at least one data symbol from each of the transmitters or may include at least one data symbol 
from each of a plurality of users. 

Other features and advantages of the present invention will become apparent from the 
following detailed description of the invention with reference to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will now be described in greater detail in the following detailed 
description with reference to the drawings in which: 



Figuro^ is a block diagram showing an example of a MTMR wireless communication 
arrangement in which plural transmitter antennas concurrently transmit respective substreams 
of symbols to plur^lreceiver antennas. 




Figure/2 is a diagram showing respective complementary cumulative distribution 




functions of the capacity C for a (1,1), (2,2) and (4,4) diversity communication system. 



Figure diagram showing respective measured bit error rates as a function of 
signal-to-noise ratio for a (2, 2) MTMR communications system using 4-PSK modulation 
with 20 Hz fading for MMSE block decision feedback sequence estimation and for ZF block 
decision feedback sequence estimation in accordance with the invention and for known 



Figureft is a diagram showing respective measured bit error rates as a function of 
signal-to-noise ratio for a (2, 2) MTMR communications system using 8-PSK modulation 
with 20 Hz fading for MMSE block decision feedback sequence estimation and for ZF block 
decision feedback sequence estimation in accordance with the invention and for known 
MMSE and ZF BLAST estimations. 

DETAILED DESCRIPTION OF THE INVENTION 

The present invention provides a MTMR architecture for detecting data transmitted 
over a wireless connection in which superimposed substreams of a signal are received and 
demodulated by concurrently estimating the substreams. The data is transmitted by multiple 
antennas or by multiple users. 

Figure 1 shows an arrangement for transmitting, receiving and detecting data 
according to the invention. 

At a transmitter location, a single, encoded data stream comprised of mN symbols is 
delivered to a demultiplexer (DEMUX) 1 02 which divides the data stream into m data 
substreams of N symbols each. The DEMUX delivers the substreams to transmitters 7x7, 




MMSE and ZF BLAST estimations. 




7x2, Txm, respectively. The transmission power is equally divided among the m transmit 
antennas. 

The transmitters 7x7 to Txm transmit the substreams over a common wireless channel, 
typically at a symbol rate 1/ T symbols/sec, with synchronized symbol timing. The 
transmitters may be quadrature amplitude modulation (QAM) transmitters, phase shift 
keying (PSK) transmitters, or other known transmitters that may be represented as a vector- 
valued transmitter where components of each transmitted m-vector are symbols dravm from a 
QAM constellation. 

Alternatively, the m data substreams are transmitted by a single transmitter antenna 
but represent m users. 

At a receiving location, a plurality of antennas 7?x7, Rx2, Rxn each receive signals 
transmitted from all of the m transmitter antennas or m users and deliver the received signals 
to a receiver 104 which concurrently demodulates all of the received signals. The n antennas 
and the receiver 104 are typically comprise a quadrature amplitude modulation (QAM) 
receiver system, though over modulation methods may be employed. 

The capacity C of the m transmitters and n receivers is defined by the following 
relation: 



C = log2 det I + — HH " bits/sec/Hz, 




where H is an w x channel matrix, /? is a spatial average of the signal-to-noise ratio (SNR) 
at each of the n receiver antenna elements. 



Figure 2 shows respective complementary cumulative distribution functions (ccdfs) of 
the capacity C for: (1) a (1,1) diversity communication system, namely a communication 
system comprised of one transmitter antenna and one receiver antenna, (2) a (2,2) diversity 
communication system, namely a communication system comprised of two transmitter 
antennas and two receiver antennas, and (3) a (4,4) diversity communication system, namely 
a communication system comprised of four transmitter antennas and four receiver antennas. 
As Figure 2 shows, the capacity increases linearly with the number of transmitter and receiver 
antennas. 

The detection operation of the receiver+t© is now described: 

A signal sample x^^ is received by an antenna j of the n receiver antennas at a discrete 

time instant k. The received signal sample is represented by a sum of the symbols 
transmitted by the m antennas according to the following relation: 

m 

where a^j are the channel impulse response coefficients of a signal path from the rth 
transmitter antenna to the yth receiver antenna, s\ is the symbol transmitted by the /th 

transmitter antenna, or from the /th user, that is received at the time instant A:, and n{ is a 

value representing additive white Gaussian noise that is received by the 7th receiver antenna 
at the time instant k. The channel impulse response coefficients a^j are estimated from 
training symbols which are embedded in the transmitted symbol sequence and whose values 
are already known to the receiver. The receiver compares the expected training symbol 



values with the values that are actually received and uses a least squared channel estimation 
to determine the values of the channel impulse response coefficient. 

The representation of a single signal sample received at the receiver antenna j may 

be extended to a stacked block data model, which represents N signal samples received at the 
receiver antenna j\ by stacking the received signals x^ into a Nx\ vector x!_ where: 

X — 5 X2 5 . . . 5 j 5 

in which T a transpose operator. The vector y may be represented as a function of the 
transmitted symbols in Nxl complex space by the relation: 

x' =A's-hn' eC'^^\ 

where is a matrix representing the channel response values of the signal paths between the 
first to wth transmitter antennas or users and the jih receiver antenna, ^ is an vector of the 
respective symbols transmitted by each of the m transmitter antennas, and nl is a vector 
representing additive white Gaussian noise that is received by the yth receiver antenna at the 
time that each of the signal samples is received. 

In a flat Rayleigh fading environment, in which there is multipath scattering with no 
delay spread, the matrix is ordinarily a block diagonal matrix that is defined in NxNm 
complex space by the relation: 




# 



namely AMs a row, Nm column matrix that is arranged as follows: 



a^j 0 0 0 
0 a^j 0 0 
0 0 a^jO 
0 0 0 «, 



• a^jO 0 0 a^jO 0 0 



0 a^j 0 0 



0 « 0 0 

mj ^ 



0 0 a^jO 0 0 a^jQ 

0 0 0 a^j 0 0 0 a^j 



The invention, however, also accounts for multipath scattering with a delay spread, 
i.e. accounts for the presence of ISI, and therefore the matrix typically includes values that 
"band" the block diagonal values. Thus, A^ is a block banded matrix with N+v rows and Nm 
columns that is arranged as follows: 



a\jQ 0 0 a\jO 0 0 aljO 0 0 

ala\,0 0 a,^.fl[',0 0 aLaLO 0 
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0 0 : - 



The above vector s is comprised of the symbols transmitted by each of the m 
transmitter antennas and is defined by the relation: 



^ = [(^')^(^0' U'Y 



in which s! is a vector representing the N symbols transmitted by the iih transmitter antennas. 
The vector s' is defined by the relation: 



The vector represents the additive white Gaussian noise and is defined by the 
relation: 



The signal samples received by each of the n receiver antennas may be then combined 
into a single vector representation as follows: 

7 = 1 



where the vector A is defined according to the relation A = ^ ^ [ A'' ] ^ [ A"' ] in which 

[A^]^ is the complex conjugate transpose of [A^]. When the above relation for x! is inserted 
into X, the following relation is obtained: 
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X 



= Xs + n € C 



Nmx\ 



where w represents the additive white Gaussian noise, ^ = X "^J^"^]^^? with a covariance 



matrix given by the relation: 



/ 1^ r 



cr'A, 



where (S' is an expectation operator, is the noise variance, and it is assumed that: 



Using the above relation for x, the values within the transmitted data symbols vector 
are determined using decision-feedback sequence estimation to satisfy a performance related 



criterion, such as a minimum mean squared error (MMSE) or a zero forcing (ZF) criterion. 
When MMSE is the performance related criterion, a minimum mean squared error-block 
linear estimator (MMSE-BLE) is derived from which a minimum mean squared error-block 
decision feedback sequence estimator (MMSE-BDFSE) is obtained. When ZF is the 
performance related criterion, a zero forcing-block linear estimator (ZF-BLE) is derived from 
which a zero forcing-block decision feedback sequence estimator (ZF-BDFSE) is obtained. 

To determine the MMSE block level estimator, the values within the matrix of the 
transmitted data symbols s are estimated by minimizing the following quadratic form Q^viSE' 



s- s 



= mm^' 



where s is an vector of estimated data symbols. 



s - s 



^~ ^) [l~^)is the magnitude 



of the vector 5-5, and s is the estimated value of s. By solving for J, the following relation is 
obtained: 



= arg min 



= (A"R:,;A+R;')-'A"Rr:i 

=(|,[a^]"[a^>.'r;')"(|,[a^]%>), 
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where is it assumed that s is distributed according to a Gaussian distribution having a zero 
mean value and a covariance matrix R^^, namely s — N (0, R^s), and n is distributed according 
to a Gaussian distribution having a zero mean value and a covariance matrix K^^, namely 

n^N(0, RJ, and R, = s{ss''} = 1. 

To determine the ZF block level estimator, the values within the vector of the 
transmitted data symbols s are estimated by minimizing the following quadratic form Qzr: 



min Qzpi^ ) = min^ \\x - A^H 

R 

||x-A. )"r-:(x-A. )}. 



= mm 



The solution for s is given by the relation: 

A 

s = ar§,minQ2f(s) 

= (a^r»->a^r;:x 

= (f[A>r[A^])"'(z[A>]%>), 
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where it is assumed that n is distributed according to a Gaussian distribution having a zero 
mean value and a covariance matrix R^n, namely n — N(0, ^nJ- 

The proof of the solutions for the MMSE block level estimator and the ZF block level 
estimator are known, such as is described in "Fundamentals of Statistical Signal Processing - 
Estimation Theory", Vol. 1, S. M. Kay, Prentice Hall, 1993. 

Improved detection performance is obtained using the MMSE block decision 
feedback sequence estimator and the ZF block decision feedback sequence estimator, instead 
of the MMSE block level estimator and the ZF block level estimator. The MMSE block 
decision feedback sequence estimator and the ZF block decision feedback sequence estimator 
are derived using a Cholesky decomposition of the inverse portions of the above expressions 
for the MMSE block decision feedback sequence estimator and the ZF block decision 
feedback sequence estimator. Namely, using the expression for the MMSE block level 
estimator, the relation 



is defined and then substituted into the above MMSE expression for i to obtain the following 
relation: 





V 7=1 



Similarly, using the ZF block level estimator, the relation 
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is defined and then substituted into the corresponding ZF expression for s to likewise obtain 
the relation: 



1 [A^] X 

\J=\ 



j 

J 



For both MMSE and ZF, the above expression for ican be muhiplied by to obtain a 
parameter z that is defined by the relation: 



( n H \ 



The solution for 5 is a continuous valued function. 

To obtain a discrete valued solution, an absolute difference value A, defined as 



A = 



J is minimized for all possible estimated values s of the transmitted symbols 



5, such as the phase shift keying (M-PSK) constellation points, to obtain the following 
expression: 
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To reduce the estimation error and increase the speed of the estimation, the discrete 
value minimization is applied to respective blocks of data symbols, rather than to the data 
symbols individually, so that an entire block of data symbols can be estimated at a time. The 
above expansion can therefore be rewritten to obtain the following cost function: 




where 



U- 



^W{i~ 1 ),W{j- 1 ) * " ^W{i-\\{ W-\)j 



I •••/ 



G C 
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S(iv-\)i 



m 

5?n N/W. 



The value W is the size of the block of data symbols, namely the number of data symbols in a 
block, and the value Q is the number of blocks in the N transmitted data symbols, namely Q = 



By substituting the above relations, the cost function can be rearranged into the 
following form: 



1." s-z 



^QQ-Q -a 



^Q-\,Q-\-Q-\ ^e-ie-e -Q-\ 



Q 


2 




Q 






••4- 


>1 



= A. 
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Then, the parameter A is minimized using the following iterative algorithm: 

(1) For the first error term A^^ in the cost function, possible values for the estimated 

symbols Sq are tested until a W sized sequence of symbols that minimizes the first 
error term A_q (sq) is obtained. This operation is referred to as sequence estimation. 

(2) Using the symbols Sq estimated in the previous step, possible values of Sq.^ are 
tested until a sized sequence of symbols is obtained that minimizes the error term 
4q-i(Sq.i, Sq). This operation is known as decision-feedback. 

(3) The above decision-feedback operation is repeated to minimize the error term 
Aq_2. Namely, using the symbols Sq and Sq., estimated in the previous step, possible 
values of Sq_2 tested to obtain a W sized sequence of symbols that minimizes the error 
term Aq.2(Sq.2, Sq.,^ Sq). 

(4) The above decision-feedback operation is repeated for each error term until the 
final error term A, is estimated, i.e, until, using the symbols S2, . . . , Sq estimated in 
the previous steps, possible values of s, are tested until a W sized sequence of 
symbols is obtained that minimizes the error term Ai(si, . . Sq). 

The number of possible symbol states that are tested for each block is Af^. Thus, the 
size of should be restricted to values that attain an acceptable detection performance while 
the complexity of processing the test values is kept within the capability of the digital signal 
processor or other processor that carries out the above iterative process. As an example, for 
eight phase state phase shift keying (8-PSK) modulation, W preferably has the value 3 or 4. 

Figure 3 shows a measured bit error rate (BER) as a function of the signal-to-noise 
ratio (SNR) for a two transmit antenna, two receive antenna, known as (2, 2), MTMR 
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communication system using 4-PSK modulation with 20 Hz fading. Five thousand 
independent time slots are averaged for each BER calculation using tw^enty five training 
symbols for each antenna used for channel estimation. The lower solid curve (identified as 
MMSE-BDFSE) represents a MMSE block decision feedback sequence estimation of the 
5 invention, and the upper solid curve (ZF-BDFSE) represents a ZF block decision feedback 
sequence estimation of the invention. The lower dashed curve (MMSE-BLAST) and the 
upper dashed curve (ZF-BLAST) represent MMSE BLAST and ZF BLAST estimations, 
respectively. The MMSE block decision feedback estimation of the invention attains a lower 
BER for a given SNR than the MMSE BLAST estimation, and the ZF block decision 
1 Ol'n feedback estimation of the invention attains a lower BER for a given SNR than the ZF 

=p BLAST estimation. 

fn 

Figure 4 illustrates the measured BER as a fiinction of the SNR for a (2, 2) MTMR 
communication system using 4-PSK modulation with 20 Hz fading. Five thousand 

ry independent time slots are averaged for each BER calculation using twenty five training 

s Li 

15Cm symbols for each antenna used for channel estimation. The lower solid curve (MMSE- 
O BDFSE) represents a MMSE block decision feedback sequence estimation of the invention, 
and the upper solid curve (ZF-BDFSE) represents a ZF block decision feedback sequence 
estimation of the invention. The lower dashed curve (MMSE-BLAST) and the upper dashed 
curve (ZF-BLAST) represent MMSE BLAST and ZF BLAST estimations, respectively. The 

20 MMSE block decision feedback estimation of the invention attains a lower BER for a given 
SNR than the MMSE BLAST estimation, and the ZF block decision feedback estimation of 
the invention attains a lower BER for a given SNR than the ZF BLAST estimation. 
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In contrast to existing techniques, such as the BLAST algorithm, the receiver scheme 
of the invention allows for the estimation of the transmitted symbols in the presence of inter- 
symbol interference (ISI). Thus, the above process may also be applied to other 
communication systems, such as those incorporating Turbo processing and orthogonal 
frequency division multiplexing (OFDM). 

In OFDM based wireless systems, for example, the same subcarriers are used in 
different transmit antennas. The invention may be used to demodulate the transmitted signal. 
The invention may also be applied to multi-user detection or multi-code detection in third 
generation, universal mobile telecommunication system (UMTS) base transceiver station 
(BTS) products. 

Although the present invention has been described in relation to particular 
embodiments thereof, many other variations and modifications and other uses may become 
apparent to those skilled in the art. It is preferred, therefore, that the present invention be 
limited not by this specific disclosure herein, but only by the appended claims. 
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